Power law eigenvalue density, scaling and critical random matrix ensembles 



K. A. MuttalilE 

Department of Physics, University of Florida, P.O. Box 118440, Gainesville, FL 32611-8440 

Mourad E.H. IsmaiQ 

Department of Mathematics, University of Central Florida, Orlando, FL 32816 

We consider a class of rotationally invariant unitary random matrix ensembles where the eigen- 
value density falls off as an inverse power law. Under a new scaling appropriate for such power 
law densities (different from the scaling required in Gaussian random matrix ensembles), we cal- 
culate exactly the two-level kernel that determines all eigenvalue correlations. We show that such 
ensembles belong to the class of critical ensembles. 
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I. INTRODUCTION 

Gaussian random matrix ensembles, introduced by 
Wigner and Dyson [1], have been studied extensively over 
half a century in the context of nuclear physics, atomic 
and molecular physics, condensed matter physics as well 
as particle physics The wide applicability results 
from the universal properties of the correlation functions 
of the eigenvalues, once they are appropriately scaled. 
Thus, e.g., the correlation functions of the eigenvalues x 
of an ensemble oi NxN hermitian matrices become, after 
proper scaling, independent of the size of the matrix or of 
the details of the microscopic distribution from which the 
matrix elements are drawn. For a Gaussian distribution, 
the density is a semicircle given by cr(x) — ^2N — x^/tt; 
this requires, for universality, the double scaling limit 



OO, 



0, u = xViV finite. 



(1) 



such that the density u(u) — constant near the origin. 
This particular double scaling is required for all Gaussian 
random matrix ensembles of different symmetry classes 
in order to have a universal large N limit, fn par- 
ticular, for Gaussian unitary ensembles, after a second 
trivial scaling Q = v^w/tt such that i7(^) = 1, it leads to 
the universal two-level kernel 
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independent of N , which gives rise to universal eigen- 
value correlations like the nearest-neighbor spacing dis- 
tribution or the number variance, commonly known as 
the Wigner distributions. 

Note that in the above double scaling limit, one is al- 
ways restricted to the eigenvalues far from the tails of the 
density. This is a highly non-trivial feature of Gaussian 
random matrix ensembles. Indeed, the eigenvalue corre- 
lations for Gaussian ensembles are deeply related to the 
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properties of Hermite polynomials 3 , and the above scal- 
ing is dictated by the asymptotic properties of Hermite 
polynomials of order N and argument x in the large N 
limit for small x. Since Hermite polynomials have a dif- 
ferent asymptotic behavior for large x in the large N limit 
[J, a different scaling is dictated for the Gaussian en- 
sembles near the edge of the semi-circle spectrum. Here 
universality is recovered after a shift x' = x — 
and a different double scaling limit near the new origin, 
namely A'^ — *■ 00, x' — > 0, with u = x'N^^^ remaining 
finite 0. In this case one gets the so called Airy kernel 
and the density decreases exponentially for large u. 

In a sense, the universality of the Gaussian random ma- 
trix ensembles is a result of an underlying central limit 
theorem, generalized to matrices Q. Indeed, the univer- 
sality requires that the eigenvalues are confined 'strongly 
enough' such that the fluctuations remain bounded. For 
example, if the probability density P{H) for an A^ x A^ 
random Hamiltonian H is given by 



P{H) cx exp[-/3TrF(iJ)], 



(3) 
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then in a more general situation where the confining po- 
tential V{H) is not a Gaussian, the density may not be 
the same near the origin. Even then, one can 'unfold' 
the spectrum by choosing a variable in which the den- 
sity is unity, and as long as the confining potential is 
'strong enough', the double scaling given by Eq. (1) al- 
ways leads to the same two-level kernel and therefore the 
same universal correlations. This is because for any such 
rotationally invariant ensemble described by a confining 
potential, the two-level kernel can be described in terms 
of polynomials orthogonal with respect to a Freud-type 
weight function e"^*-^-*, with polynomial V, (just as the 
Hermite polynomials arise for Gaussian ensembles where 
the potential is V{x) = Again, the eigenvalue 

correlations are deeply related to the properties of those 
orthogonal polynomials. However, all Freud-type orthog- 
onal polynomials have qualitatively similar asymptotic 
behavior in the large A^ limit. Therefore the same dou- 
ble scaling limit leads to universality for all such differ- 
ent confining potentials. Conversely, an arbitrary scal- 
ing (with e.g. an arbitrary power law for A^) does not 



2 



lead to a universal iV-independent kernel for rotationally 
invariant random matrix ensembles associated with the 
classical orthogonal polynomials. 

On the other hand, it is known that the universality 
does break down if the confinement potential V{H) grows 
extremely slowly with H [7|, namely 



V{H) 



1 
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In^ H, iJ > 1. 



(4) 



This leads to the so called 'critical' ensembles [8[, charac- 
terized by a one-parameter generalization of the Gaussian 
two- level kernel, given by J, 9] 
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The eigenvalue correlations for such systems for various 
values of the parameter 7 have been studied in the con- 
text of the Anderson transition in disordered quantum 
conductors [IQ,] • The density of the eigenvalues with the 
weak confinement potential (4) turns out to be a con- 
stant independent of N in the N 00 limit, and the 
scaled variables are simply C, — x/j. Note that the criti- 
cal ensembles tend to the Gaussian ensembles in the limit 
7 — > 0, in which case one is again restricted to the corre- 
lations of the eigenvalues near the origin. However, the 
scaling does not involve N. This difference in the scaling 
behavior is the reason why critical ensembles could not 
be studied using the classical orthogonal polynomials. 

A second case where a different scaling beyond the 
scope of the Freud-type orthogonal polynomials is re- 
quired is when the density of the eigenvalues x falls off as 
a power law for large x. This indeed happens for a large 
class of random matrix ensembles of recent interest, rel- 
evant for e.g. finance or scale free systems that requires 
considering power law distributions with fat tails pl Ul3 | . 
In a pioneering study, Cizeau and Bouchaud (CB) Q in- 
troduced the so called 'Levy matrix ensembles' where the 
matrix elements are drawn from a power law distribution 
according to P{Hij) with 



Hi, > 1. 



(6) 



Here Hq is of order Hij and the parameter /i > 0. 
The eigenvalue density in this case was shown to fall off 
with the same exponent, as l/x^~^^. In particular for 
< fi < 1, the distribution has diverging variance and 
hence the underlying central limit theorem must now be 
modified according to the theorems of Levy and Gne- 
denko (l3| . Indeed, several numerical work have shown 
(T3I [1^ that the universality of the Gaussian ensembles 
does break down when the density of eigenvalues follows 
a power law with exponent less than 2, the signature of 
which is apparent in the numerically obtained nearest- 
neighbor spacing distribution and the number variance. 
The eigenvalue densities for certain Levy ensembles have 
been obtained analytically [12], but so far it has not been 



possible to evaluate the two-level kernel which gives the 
two or higher level eigenvalue correlations. Since the uni- 
versality requires 'unfolding' where the density is made 
uniform, it is obviously of more interest to obtain the two- 
level kernel for a Levy-like ensemble, specially where the 
microscopic distribution leads to a power law eigenvalue- 
density which falls off with an exponent less than 2. 

Clearly, if the density of eigenvalues falls off with a 
power law and we are interested in the properties of the 
large eigenvalues in the tails, any possible universality 
of the associated random matrix ensembles cannot be 
expected in the same double scaling limit used for the 
Gaussian ensembles, either at the origin or near the semi- 
circle edge. In @, CB proposed a double scaling limit 
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x/N" finite. 
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where a > is related to the power law exponent, in 
which the Levy ensembles might become universal. We 
will call this the CB scaling. As mentioned above, this 
scaling cannot be obtained naturally within the scheme 
of classical or Freud- type orthogonal polynomials. This 
means that the most powerful approach to study eigen- 
value correlations of random matrix ensembles, namely 
the method of orthogonal polynomials [H, seem to be 
inapplicable for Levy-type random matrices. It is for 
this reason that progress has been slow in detailed ana- 
lytic studies of eigenvalue correlations of random matri- 
ces with power law densities. 

On the other hand, although in most studies it has 
been implicitly assumed that a power law density of 
eigenvalues requires choosing matrix elements from a 
power law distribution of the type (6), we find that at 
least for the inverse power law with exponent 1, this is 
not the case. Indeed, we find that a weakly confined log- 
squared potential that gives rise to the critical ensemble 
with a constant density for an A'^-independent scaling also 
gives rise to a inverse power law density under the CB 
scaling. It turns out that the weakly confined log-squared 
potential is exactly solvable even in the new CB scaling 
regime and the two-level kernel in this scaling regime has 
a well-defined limit as A^ ^ cx3. 



II. THE MODEL 

In the present work we use the solvability of the log 
squared potential to obtain analytically, exactly and ex- 
plicitly, the two-level kernel for a random matrix ensem- 
ble with density a{x) oc 1/x in the CB double scaling 
limit given in Eq. (7) where a is an arbitrary positive 
parameter. Specifically, we use the model characterized 
by the confinement potential 
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where g is a parameter, < g < 1. The joint probability 
distribution of the N eigenvalues can then be written as 
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P{{x}) ^ CN{q)l[{x. 

i<j 



4=1 



-Vix\q} 



(9) 



where Ciy{q) is a constant and the product Y[i<j(^i~^j)'^ 
is the standard Vandermonde factor for unitary ensem- 
bles. The orthogonal polynomials corresponding to the 
weight function 



WH{x\q) = q^/^ 



TT In q 



'V{x\q) 



(10) 



for the above confining potential, in terms of which the 
two-level kernel can be written down exactly, are the 
Ismail-Masson q^^-Hermite polynomials 15] first consid- 
ered in the context of random matrices in ^ , where only 
an N-independent scaling was considered. The reason 
for the choice of the potential ([5]) is that a very general 
asymptotic relation for the corresponding Ismail-Masson 
polynomials has been obtained recently p^ . such that 
the new CB scaling can be implemented explicitly with- 
out compromising the solvability of the model. In par- 
ticular, for 



Xn{t,u) = -{q 



q'^'/u) 



(11) 



and < i < 1/2 [Tj] it has been shown that the normal- 
ized Ismail-Masson polynomials hn{xn\q) has the large 
n limit given by 

wh (sinh hiu\q)u^'^ 



WH{Xn) hn{Xn) 



(-l)L™/2J(5;^)^y(^ 

q^^^U^'ei-u'q^'^^q) (12) 



Here {q; q)n = i^ — q^), is an integer and < A < 1 
has been defined by the relation m + \ = n{\ — 2t). We 
have also defined 2a — x(™) + A where x(™) is 0(1) for 
m =even(odd). The function 



e(z,g)= £ q^' 



is defined for z 7^ for any complex z. 
Let us define 



q -~ 



7 > 0. 



(13) 
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Note that for ^ 00 and large u, we can choose the 
parameters 7 and t to satisfy the CB scaling for arbitrary 
a. The choice 



-fNt = a In iV, 



(15) 



when used in Eq. (11), gives precisely the CB scaling u = 
x/N°'. Note that while for finite < i < 1/2 the above 
condition requires 7 = 2aln(A^)/A^ or equivalently 
q ^ 1, we are free to choose t ~ ln(A^)/A^ such that 
the scaling relation can be satisfied for finite 7. In the 
following, we will outline the derivation of the two-level 
kernel for 7 <C tt^ . 



III. THE TWO-LEVEL KERNEL 

The two-level Ismail-Masson kernel defined in terms 
of the orthonormal polynomials hn{x\q) is given by 

KN{x,y) \/'^H{xN)wH{yN)J2kJo^ hk{x\q)hk{y\q). 

For large N it can be written, using Christoffel-Darbeau 
formula, as 



kN 



KN{x,y) = \/wh{xn)wh[vn)-, x 

hN{xN\q)hN^i{yN\q) - hN{yN\q)hN-i{xN\q) 

x-y 



(16) 



where kjs; is the coefficient of the term x'^ in the polyno- 
mial of order N . In Eq. (16) we can replace by TV — 1 
to obtain the asymptotic expression for the normalized 
polynomials hN^i{xN~i\q)- However, the Christoffel- 
Darbeaux formula involves hi^-i{xiq\q)- Similarly, we 
need expression for the weight factor WH{xN-i\q) in 
terms of WH(xN\q)- We obtain these results by exploiting 
the scaling relation 



XN-i{t,q *u) = XNit,u), 



(17) 



which follows from Eq. (11). This gives, in terms of the 
scaled variables v = y/N"-, 



V^HiVN) hN-i{yN 



^ Wh (sinh hi{q~'^v)\q) 



(-l)L'"'/2j(g;q)^^(q;g)^_, 

((?-*«)2^+2*g(^-i)/4q(6+*)^e(-«2q2f ^ q) (18) 

where m' + X' = {N — 1)(1 — 2t) and we have defined 
26 = x(m') + X' - 2t. 

The functions Q{z,q) appearing in the above equations 
cannot be evaluated easily. However, for q near 1 (or 
7 <C TT^) which is the region we would be interested in 
(this allows us to go to the Gaussian limit q 1), we 
can use the imaginary transformation [181] 



/ TT In^ w TV In iv^ 

Q{w,q) — . —e '^-< 8(e ^1 ,p); p — e -< . (19) 
V 7 



In this case the series in powers of p = e~'^^^^ converges 
rapidly for 7 <C tt^ and can be approximated quite well 
by keeping only the dominant term. The transformation 
allows us to rewrite the kernel as 

(x - y) K{x, y) = A(q)e'^ ^ pkl+kl-k^-k, ^ 

ki ,k2 

[pzif' {pz2f'' [1 - e-^2^("-'')('=i-'=^)] (20) 

where pzi = e'^'^I"'^ , pz2 = e^^'^I^"^ and A{q) is 
a constant independent of We now make our first 
approximation, namely that we keep only the two terms 
ki — Q,k2 — 1 and A;i = 1, A:2 = in the double sum in 
Eq. (20). This neglects terms of order or higher powers 
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of p. Using the fact that the product (-1)L™/2J + L'n72j 
sm[7r(a — b)] — 1 for all a, 5, we finally obtam 



K{x, y) = Koq 



sin[^ In 2il 

nt 1-7 



(21) 



Thus it shows that rotationally invariant random matrix 
ensembles with inverse power law density a{x) oc 1/x 
belong to the class of critical ensembles, albeit after a 
logarithmic variable transformation. 



where 



r2 27r ,/„ e 



(22) 



In order to estimate Kq, we use the Euler identity [T3| 

(9;9)oo = nr=i(i-9") = Er=-oo(-i)"'?"^'"+'^/'- The 

sum can be evaluated using the Poisson summation for- 
mula to yield 



A;— — cx 



e-lT(2fe+i) -^(2fe+i), (23) 



Again we use the fact that 7 ^ tt^ and approximate the 

series by keeping only the fc = and A: = — 1 terms. 
— ^2 

yq"^/^"'e^"3^. This then leads to 

Kn = 1/tt. 

IV. RESULTS AND DISCUSSION 

In terms of the scaled variables u = x/N", v = y/N", 
the two-level kernel in the large N limit is given by 

1 sin[- In -1 

i^^(u,w;7) = ^ 0<7<^' (24) 



with density 



TT U — V 



a(u) = {u, u) 



7u 



(25) 



Thus the kernel is not universal, and depends on one 
parameter that cannot be scaled out. Surprisingly, once 
unfolded by changing variables to 



C = — Inu, ?7 = — Inw 

7 7 



(26) 



such that the density becomes unity, a{(^) = 1, the scaled 
kernel becomes identical to the critical kernel of Eq. (5), 



(27) 



It is important to emphasize that although the current 
model belongs to the same class of critical ensembles, it 
differs from the A^-independent scaling model in two im- 
portant ways. First, the A^-dependent scaling introduces 
non-trivial A^-dependent logarithmic potential in addi- 
tion to the log-squared potential of the A^-independent 
scaling model. Second, the logarithmic variable transfor- 
mation changes the interaction between the eigenvalues. 
It is therefore not obvious that the present model should 
lead to a critical ensemble. In the variable where both 
models have weakly confining log-squared potential, the 
kernels for the CB scaling and the A^-independent scaling 
are given by Eqs (24) and (5), respectively. 

It has been conjectured before that an inverse 
power law density in a rotationally invariant random ma- 
trix ensemble approaches the same behavior as that for 
a weakly confined log-squared potential considered here. 
In this work we show explicitly by calculating exactly 
the two-level kernel that the inverse power case is a crit- 
ical ensemble under the CB scaling. However, it differs 
from the case of the weakly confined log-squared poten- 
tial with A^-independent scaling in important ways; in 
the variable where the present kernel becomes identical 
to that for A^-independent scaling, the confining potential 
is in fact Gaussian. 



The Ismail-Masson polynomials go back to the classi- 
cal orthogonal polynomials in the limit q ^ 1. One can 
therefore hope to implement the scaling for the Ismail- 
Masson polynomials and then take the g — s- 1 limit to 
obtain the corresponding scaling for the Gaussian en- 
sembles. It turns out that the inverse power density 
a{u) cx 1/u is very special and exists only for q < 1. 
The unfolded kernel with constant density of course has 
a well defined g — > 1 limit which is just the sine kernel of 
the Gaussian ensemble. It should be possible to obtain 
the two-level kernel for arbitrary density power law with 
CB type scaling within our scheme, but it is not clear if a 
well-defined N ^ 00 limit exists for the resulting kernel. 
Here we presented only the inverse power density case, 
which we show to be critical. 



[1] E.P. Wigner, Proc. Cambridge Philos. Soc. 47, 790 

(1951); F.J. Dyson, J. Math. Phys. 3, 140 (1962). 
[2] For recent reviews and applications see T. Guhr, A. 



Mueller-Groeling and H. Weidenmueller, Phys. Rep. 299, 
189 (1998); Special issue: Random Matrix Theory, J. 
Phys. A 36, (2003), Guest Eds. P.J. Forrester, N.C. 



5 



Snaith and J.J.M. Verbaarschot. 
[3] M.L. Mehta, Random Matrices, 2nd Ed, Academic Press, 
1991. 

[4] G. Szego, Orthogonal Polynomials, 4th Ed., AMS, RI, 
1975. 

[5] P.J. Forrester, Nucl. Phys. B 402 [FS], 709 (1993); C.A. 
Tracy and H. Widom, Commun. Math. Phys. 159, 151 
(1994). 

[6] P. Cizeau and J. P. Bouchaud, Phys. Rev. E 50, 1810 
(1994). 

[7] K.A. Muttahb, Y. Chen, M.E.H. Ismail and V.N. 

Nicopoulos, Phys. Rev. Lett. 71, 471 (1993); C. Blecken, 

Y. Chen and K.A. Muttalib, J. Phys. A 27, L563 (1994). 

For a review, see K. A. MuttaUb, Y. Chen and M. E. H. 

Ismail, in Symbolic computation, number theory, special 

functions, physics and combinatorics Eds. F. Garvan and 

M. Ismail, Kluwer Academic (2001). 
[8] V.E. Kravtsov and K.A. Muttalib, Phys. Rev. Lett. 79, 

1913 (1997); E. Bogomolny, O. Bohigas and M.P. Pato, 

Phys. Rev. E 55, 6707 (1997); E.B. Bogomolny, U. Ger- 

land and C. Schmit, Phys. Rev. E 59, R1315 (1999); 

A.M. Garcia-Garcia and J.J.M. Verbaarschot, Phys. Rev. 

E 67, 046104 (2003). 
[9] M. Moshe, H. Neuberger and B. Shapiro, Phys. Rev. Lett. 

73, 1497 (1994); 
[10] A.D. Mirlin, Y.V. Fyodorov, F.M. Dittes, J. Quezada 

and T.H. Seligman, Phys. Rev E 54, 3221 (1996); S.M. 



Nishigaki, Phys. Rev. E, 59, 2853 (1999); A.D. Mirlin 
and F. Evers, Phys. Rev. B 62, 7920 (2000). 

[11] L. Laloux, P. Cizeau, J. P. Bouchaud and M. Potters, 
Phys. Rev. Lett. 83, 1467 (1999); A. Moustakas et al. 
Science 287, 287 (2000); M. Araujo, E. Medina and E. 
Aponte, Phys. Rev. E 60, 3580 (1999). 

[12] Z. Burda et al, Phys. Rev. E 65, 021106 (2002); 

[13] J.P. Bouchaud and A. Georges, Phys. Rep. 195, 127 
(1990). 

[14] J. Choi, unpublished. 

[15] M.E.H. Ismail, Classical and Quantum Orthogonal Poly- 
nomials in one Variable, Cambridge University Press, 
Cambridge, 2005. 

[16] M.E.H. Ismail and R. Zhang, IMRN 2006, ID 83274, 1 
(2006). 

[17] The asymptotic relation for the Ismail-Masson polyno- 
mials has been obtained in [l^ for all t > 0. Here we use 
only the result for t < 1/2. 

[18] E.T Whittaker and G.N. Watson, A Course of Modem 
Analysis, 4th Ed., Cambridge University Press, 1927. 

[19] See e.g. G. Gasper and M. Rahman, Basic Hypergeo- 
metric Series, Cambridge University Press, Cambridge, 
1990. 

[20] A.C. Bertuola, O. Bohigas and M.P. Pato, Phys. Rev. E 
70, 065102(R) (2004). 



